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Does the Sun Appear Brighter at Night in Neutrinos? 
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We calculate accurately the number of solar neutrino events expected as 
a function of solar zenith angle, with and without neutrino oscillations, for 
detectors at the locations of Super-Kamiokande, SNO, and the Gran Sasso 
National Laboratory. Using different earth models to estimate geophysical 
uncertainties, and different solar models to estimate solar uncertainties, we 
evaluate distortions predicted by the MSW effect in the zenith angle distribu¬ 
tions of solar neutrino events. The distortions are caused by oscillations and 
hy u — e interactions in the earth that regenerate Ue from or v-,-. We show 
that the first two moments of the zenith-angle distribution are more sensitive 
to the small mixing angle MSW solution than the conventionally studied day- 
night asymmetry. We present iso-sigma contours that illustrate the potential 
of Super-Kamiokande, SNO, BOREXINO, ICARUS and HERON/HELLAZ 
for detecting the earth regeneration effect at their actual locations (and at 
the equator). MSW solutions favored by the four pioneering solar neutrino 
experiments predict characteristic distortions for Super-Kamiokande, SNO, 
BOREXINO, and ICARUS that range from being unmeasurably small to 
> 5(7 (stat) after only a few years of observations. 
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I. INTRODUCTION 


Four operating solar neutrino experiments (Chlorine [Q, Kamiokande [Q|, GALLEX [Q) 
and SAGE [Q have detected neutrinos from nuclear fusion in the interior of the sun with ap¬ 
proximately the numbers and energies expected from standard solar models . Moreover, 
sound speeds calculated from the standard solar models agree with the helioseismologically 
determined sound speeds to a rms accuracy of better than 0.2% throughout essentially the 
entire sun [^. 

Nevertheless, quantitative discrepancies have persisted for almost three decades between 
the predictions of the standard solar models and the observations of solar neutrino exper¬ 
iments Ml- Several suggested modihcations of neutrino properties provide excellent hts 


to the existing solar neutrino data [11 


Are there potential “smoking gun” indications of new physics? Yes, the most popular 
neutrino physics solution, the Mikheyev-Smirnov-Wolfenstein (MSW) effect |^, predicts 
several characteristic and unique phenomena. The MSW effect explains solar neutrino ob¬ 
servations as the result of conversions in the solar interior of z/g produced in nuclear reactions 
to the more difficult to detect or z/,-. 

Potentially decisive signatures of new physics that are suggested by the MSW effect 
include observing that the sun is brighter in neutrinos at night (the ‘earth regeneration 
effect’') 


], detecting distortions in the incident solar neutrino energy spectrum |]Ig, and 
observing that the flux of all types of neutrinos exceeds the flux of just electron type neutrinos 
MM- A demonstration that any of these phenomena exists would provide evidence for 
physics beyond the minimal standard electroweak model. 

The regeneration effect is an especially powerful diagnostic of new physics since no differ¬ 
ence is predicted between the counting rates observed during the day and at night (or, more 
generally, any dependence of the counting rate on the solar zeith angle) by such popular al¬ 
ternatives to the MSW effect as vacuum oscillations [|T^], magnetic moment transitions ||20|| , 
or violations of the equivalence principle 
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In this paper, we investigate the sensitivity of new solar neutrino experiments, Super- 
Kamiokande [0, SNO [g|], ICARUS 0, BOREXINO [g|], HERON [|§ and HELLAZ 
to the earth regeneration effect. The MSW effect predicts that, for certain values of the 


neutrino masses and mixing angles, v — e interactions in the earth (at night) may convert 

or Uj. from the sun back into the more easily detectable Vf.. 

An accurate evaluation of the systematic signihcance of experimental results will require 
the detailed Monte Carlo simulations that will be carried out by the experimental collab¬ 
orations; the collaborations will determine the best estimates and uncertainties for all the 
quantities that affect the experimental result. These results will then be analyzed using 
computer codes that include the experimental details and which make use of optimal sta¬ 
tistical techniques such as maximum likelihood analysis. In the absence of detailed Monte 
Carlo simulations of the experimental characteristics (yet to be determined), we estimate 
in this paper the statistical signihcance of expected results by comparing the predictions of 
various MSW scenarios with respect to the no-oscillation scenario including only statistical 
errors and analyzing the results with a statistic. 

The reader who wants to see the approximate power of the new experiments can turn 
immediately to Fig. ^ which shows the signihcance level (statistical errors only) with which 
new solar experiments could detect the regeneration ehect. Experience with the operating 
experiments over a period of years may be necessary to determine the size of the systematic 
errors. 

This paper is organized as follows. In Sec. |1| we summarize the general features of z/g 


regeneration in the earth. We then describe in Sec. |T| the diherent models of the earth 
used for estimating the uncertainties in the numerical calculations due to uncertainties in 
the density prohle and the chemical composition of the earth’s interior. In the detailed 
calculations that follow, we use the differences between the results obtained from the different 
models of the earth to determine the geophysical uncertainties in the MSW predictions. 

After these preliminary considerations, we determine in Sec. E3 the regions in the mixing 
angle and mass difference plane, sin^ 29 - Am^, that are allowed by the latest solar neutrino 
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data, taking into acconnt the influence of earth regeneration on the predicted counting 
rates in the chlorine, Kamiokande, GALLEX, and SAGE experiments. In addition to the 
familiar large mixing angle and small mixing angle MSW solutions, we And an additional 
LOW solution (large mixing angle, smaller neutrino mass difference) that has an acceptable 
confidence level only when earth regeneration is included in the calculations. 

We then investigate the sensitivity of next-generation solar neutrino experiments to the 
earth regeneration effect. We begin by defining and calculating the zenith-angle exposure 
function in Sec. 0. This function depends only on the location of a neutrino detector on 
the surface of the earth; it is independent of the characteristics of the detector. We also 
calculate in Sec. 0 the distorted zenith angle distribution that takes account of regeneration 
in the earth. 

The results for MSW regeneration given previously in the literature involve making 
approximations either in the model description of the earth or in the calculation of the 
average survival probability after regeneration, or both. Instead, we integrate numerically 
the differential equations describing the evolution of the neutrino states in traversing an 
accurate model of the earth, thereby avoiding the necessity of arguing that an approximation 
scheme is sufficiently accurate. In several tables in this paper, we present numerical results to 
a precision of 0.01%, an accuracy much higher than can be measured experimentally. These 
precise numerical predictions are given in order to illustrate the small effect on measurable 
quantities of some of the systematic differences. 

We introduce in Sec. 0 the first two moments of the zenith-angle distribution of neu¬ 
trino events and calculate the dependence of the moments on neutrino parameters for 
the new solar neutrino experiments: Super-Kamiokande, SNO, IGARUS, BOREXINO and 
HERON/HELLAZ. For comparison, we calculate accurately in Sec. ivg the conventional 
day-night asymmetry; we present values of the day-night asymmetry for the new solar neu¬ 
trino detectors mentioned above. 

Which characterization is more sensitive to new physics: the moments of the zenith-angle 
distribution or the day-night asymmetry? We show in Sec. [VIII| that the moments are more 
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sensitive to the small mixing angle MSW solntion and the day-night asymmetry is more 
sensitive to the large mixing angle solntion. Althongh less statistically powerfnl than the 
full Monte Carlo simnlations that will be carried ont by the experimental collaborations, 
the analyses nsing, e. .g, moments or the day-night asymmetry can be carried ont qnickly 
by theoreticians interested in determining whether new particle physics scenarios can be 
tested by the experiments or whether they are already inconsistent with data that have 
been pnblished. For completeness, we present in Sec. the moments of the electron recoil 
energy spectrnm for Snper-Kamiokande and SNO that were compnted inclnding the earth 
regeneration effect. 

We then show in Sec. ^ that the MSW predictions for regeneration in the earth depend 
only slightly on the adopted density prohle of the earth, the chemical composition of the 
earth, and the details of the solar model. 

Following the snggestion of Gelb, Kwong, and Rosen |^, we calcnlate in Sec. ^ the 


increase in the sensitivity to the regeneration effect that conld be achieved by bnilding 
detectors at the eqnator. We discnss and snmmarize onr main resnlts in Sec. ra¬ 
il. THE EARTH REGENERATION EFFECT 


We work in a two-nentrino mixing scheme involving v^. (prodnced in the snn) and either 
or Vt (prodnced by oscillations). Soon after Mikheyev and Smirnov snggested the MSW 


effect as a possible solntion of the solar nentrino problems, several anthors pointed ont 
in that day-night variations of the event rates in solar nentrino detectors conld provide 
spectacnlar conhrmation of the MSW effect and thns of new physics. 

The MSW solntion of the solar nentrino problems reqnires that electron nentrinos pro¬ 
dnced in nnclear reactions in the center of the snn are converted to mnon or tan nentrinos 
by interactions with solar electrons on their way from the interior of the snn to the detector 
on earth. The conversion in the snn is primarily a resonance phenomenon, which occnrs at a 
specihc density that corresponds to a dehnite nentrino energy (for a specihed nentrino mass 
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difference). 

During day-time, the higher-energy neutrinos arriving at earth are mostly (or 
with some admixture of z/g. At night-time, neutrinos must pass through the earth in order 
to reach the detector. As a result of traversing the earth, the fraction of the more easily 
detected z/g increases because of the conversion of z/^ (or p^-) to z/g by neutrino oscillations. 
For the small mixing angle MSW solution, interactions with electrons in the earth increase 
the effective mixing angle and enhance the conversion process. For the large mixing angle 
MSW solution, the conversion of z/^ (or Ut) to z/g occurs by oscillations that are only slightly 
enhanced over vacuum mixing. This process of increasing in the earth the fraction of the 
neutrinos that are z/g is called the “regeneration effect” and has the opposite effect to the 
conversion of z/g to z/^ (or u^-) in the sun. 

Because of the change of neutrino type in the earth, the MSW mechanism predicts that 
solar neutrino detectors should generally measure higher event rates at night than during 
day-time. 

Figure illustrates a schematic view of a solar neutrino detector at the geographic 
latitude, 0. Since the earth is spherically symmetric to 0(10“^'®), it is sufficient to consider 
the cross-section slice shown in the £gure.[| Two lines determine the geometry: one line 
defines the zenith direction, and the other line is the trajectory of the neutrino. The zenith 
angle a (0° < a < 180*^) between these two lines specihes the neutrino trajectory in the 
earth. The survival probability depends on the neutrino oscillation parameters, Am^ and 
sin^ 20, on the neutrino energy, E, and on the path (i.e., a) the neutrino travels through 
the earth. Since a changes due to the apparent motion of the sun, the neutrino survival 
probability should change with time as well, resulting in an asymmetric distortion of the 
angular distribution of events. 


*The polar-equatorial asymmetry of the earth changes the predicted MSW regeneration effects 
by less than or of order 0.3%, which is not observable with existing or planned detectors. 
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Real-time detectors, which record the times at which neutrinos interact within the detec¬ 
tor, are best suited for studying the earth regeneration effect. In radiochemical experiments 
the time of detection is poorly known, since a typical run usually lasts between several weeks 
(gallium) and several months (chlorine). [| 

Kamiokande, a real-time neutrino electron scattering detector, did not see any signal for 
the earth regeneration effect. The Kamiokande collaboration used this non-observation to 
rule out an important region in parameter space for which the predicted day-night asymme¬ 
try, or zenith-angle dependence, is large [^. However, the sensitivity of the Kamiokande 
detector was insufficient to probe the full Am^ - sin^ 29 parameter space for which there 
might be an appreciable day-night effect, measurable by Super-Kamiokande or SNO. 

There are several calculations in the literature pIH35[| of the expected magnitude of the 
regeneration effect in future experiments. Different groups of researchers have used different 
models of the earth in their calculations. No quantitative estimate has been made previously 
of the sensitivity of the measurable quantities to the adopted density prohle of the earth; 
each group has typically presented results using a specihc density prohle, often not the best 
available prohle. In the subsequent sections, we describe direct numerical computations of 
the earth regeneration ehect for six diherent models of the earth. Thus we quantify the 
dependence of the calculated characteristics of the earth regeneration ehect on the model 
of the earth and exhibit the corresponding uncertainties which turn out to be rather small. 


The density prohles in the six models of the earth are described in Sec. 111. 


III. EARTH MODELS 

The MSW ehect in the earth depends upon the electron number density as a function 
of radius. In this section, we determine best-estimates and a range of uncertainties for the 


iThe possibility of making extractions twice a day and counting separately day-time and night¬ 
time samples in the chlorine and iodine experiments has been considered carefully. 
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total mass density and for the chemical composition. We use the best available earth models 
and chemical composition for most of the calculations performed in this paper, but we also 
carry out calculations for hve older mass models in order to determine the uncertainties in 
the predicted MSW effects that arise from uncertainties in the model of the earth’s density 
prohle. We use a best-estimate chemical composition for the core that is inferred from the 
seismological measurements. To test the sensitivity of the MSW predictions to the assumed 
chemical composition of the earth, we make extreme assumptions that maximize or minimize 
the average charge to mass ratio and carry out calculations also for these extreme cases. 


A. Density Profiles 


The density distribution inside the earth is known with a precision of a few percent 
3^ . A large set of seismic measurements has been used to obtain the most accurate model, 
PREM (the Preliminary Reference Earth Model), for the earth’s density distribution. 


We will use the PREM model for all of our best-estimate calculations. This model has also 
been used by Lisi and Montanino as the basis of their recent analytic study of earth 
regeneration. Other models are described in Refs. p8| and . To determine the sensitivity 
to the assumed density prohle, we have performed calculations with a representative set of 
six different earth models, all spherically symmetric and with the same radius, = 6371 
km. 

Figure shows the density prohles of the six earth models. The density distributions 
in these models are divided into hve zones: a) a crust with a thickness of a few tens of 
kilometers, b) an upper mantle extending down from the crust to about 1000 km, c) a lower 
mantle down to about 2900 km, d) an outer core between 1250 km and 3480 km from the 
center, and e) an inner core of radius ~ 1220 km. The density changes abruptly between 
the inner and outer core, and also at the border between the lower mantle and the outer 
core. The positions of these abrupt changes are known with an accuracy of better than a 
percent from seismological data. 
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Table | compares the mass and moment of inertia that we have compnted for each of 
the six earth models with the measured values. The recent models given in and 
reproduce the total mass and moment of inertia with excellent precision. The older models 
P0|^3[ give slightly worse hts to the mass and moment of inertia. The models listed in the 
third to sixth row in Table | are in conflict with seismological measurements. 

The set of six models represents a sample that allows for variations of the density distri¬ 
bution larger than the uncertainties in the PREM model. As shown in Sec. ^ (and the last 
four columns of Table |), the large differences between the six density models produce rela¬ 
tively small (but not always negligible, see Sec. 0) changes in the predicted characteristics 
of the earth regeneration effect for the Super-Kamiokande, SNO, BOREXINO, ICARUS, 
and HERON/HELLAZ experiments. 

Figure ^ shows that the largest differences between the six models of the earth are in the 
core, below 2000 km. The operating solar neutrino experiments, and those currently under 
development, are located at relatively high northern latitudes. Solar neutrinos that are 
detected in these experiments never cross the inner core of the earth. Among the real-time 
experiments that are currently operating or are under construction, Super-Kamiokande is 
most sensitive to the core distribution. Nevertheless, the fraction of a year during which the 
neutrinos cross the outer core at the Kamioka site is small (~ 7%). 


B. Chemical Compositions 


Measurements of the propagation of seismological waves in the earth’s interior and studies 
of the properties of minerals under high pressure, have been combined to determine the 
chemical composition of the earth’s interior with relatively high accuracy Using the 

results of references we adopt a best-estimate charge to mass ratio, Z/A, of 0.468 for 

the core (83% Fe, 9% Ni, and 8% light elements with Z/A = 0.5) and 0.497 for the mantle 
(41.2% SiOs, 52.7% MgO and 6.1% FeO). 

A lower limit for the charge to mass ratio in the core is 0.465, which corresponds to 
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assuming a composition of 100% iron. From the seismic and mineral data, geophysicists 
have concluded that the minimum amount of iron in the core is 80%. We determine 


a maximum value of Z/A = 0.472 in the core by assuming a composition of 80% iron and 
20% light elements. The total range of the electron number density due to the imperfectly 
known composition of the core is about 1.5%. 

The chemical composition in the mantle is believed known to about 1% (see ref. [@). 
We consider here variations of —1% and —2%. The value of Z/X in the mantle cannot 
be increased signihcantly above the standard value of 0.496 because that would require the 
presence of a large amount of hydrogen in the mantle. 


IV. AVERAGE EVENT RATES AND MSW SOLUTIONS 


Including the earth regeneration effect, we have calculated the expected one-year aver¬ 
age event rates as functions of the neutrino oscillation parameters Xrn? and sin^ 26 for all 
four operating experiments which have published results from their measurements of solar 
neutrino event rates. These include the chlorine experiment, Kamiokande, GALLEX and 


SAGE. Thus we update our previous results given in Ref. , in which the earth effect was 
neglected. We take into account, as before, the known threshold and cross-section for each 
detector. In the case of Kamiokande, we also take into account the known energy resolution 


20 % la at electron energy 10 MeV) and trigger efficiency function 


We hrst calculate the one year average survival probability, Pse, for a large number of 
values of Xw? and sin^ 26 using the method described in Appendix A. Then we compute 
the corresponding one year average event rates in each detector. We perform a analysis 
taking into account theoretical uncertainties and experimental errors as described in . 


Table y summarizes the reported mean event rates from each detector. We obtain 
allowed regions in Xrn? - sin^ 26 parameter space by Ending the minimum and plotting 
contours of constant = xLin + where = 5.99 for 95% G.L. and 9.21 for 99% 
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C.L.|. 

The best fit is obtained for the small mixing angle (SMA) solution: 


Am^ = 5.0 X 10■%V^ 

(la) 

sm‘^29 = 8.7 x 10“^ 

(lb) 

which has a Xmin — 0.25. There are two more local minima of y^. 

known large mixing angle (LMA) solution occurs at 

The best £t for the well 

= 1.3 X 10“%V^ 

(2a) 

sin^ 29 = 0.63, 

(2b) 


with Xmin = There is also a less probable solution, which we refer to as the LOW 
solution (low probability, low mass), at 


Am^ = 1.1 X 10 ^eV^, 

(3a) 

sm‘^29 = 0.83. 

(3b) 


with Xmin = 6.9. The LOW solution is acceptable only at 96.5% C.L. 

Figure ^ shows the allowed regions in the plane dehned by Am? and sin^ 29. The C.L. is 
95% for the allowed regions of the SMA and LMA solutions and 99% for the LOW solution. 
The black dots within each allowed region indicate the position of the local best-fit point 
in parameter space. The results shown in Fig. ^ were calculated using the predictions of 
the 1995 standard solar model of Bahcall and Pinsonneault §], which includes helium and 
heavy element diffusion; the shape of the allowed contours depends only slightly upon the 


assumed solar model (see Fig. 1 of ref. |TI 


The results given here differ somewhat from those given earlier in ref. , both because 
we are now including the regeneration effect and also because we are using more recent 
experimental data for the pioneering solar neutrino experiments. Comparing the results 


iThe C.L. in this paper are always for two degrees of freedom. The values of are not reduced. 
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given in Eqs. ®-0 and Figs. ^ with the corresponding allowed regions obtained for the 
same input neutrino experimental data but without including the earth effect shows that 
terrestrial regeneration changes only slightly the best-£t solutions for the SMA solution 

5% in Am^ and sin^26*) and the LMA solution (< 10% in Am^ and sin^26'). The 
values of Xmin signihcantly affected. The LOW solution is acceptable only if the 

regeneration effect is included; otherwise, the LOW solution is ruled out at 99.9% C.L. 

Figure § compares the computed survival probabilities for the day (no regeneration), the 
night (with regeneration), and the annual average. These results are useful in understanding 
the day-night shifts in the energy spectrum that are computed and discussed in Sec. The 
results in the hgure refer to a detector at the location of Super-Kamiokande, but the differ¬ 
ences are very small between the survival probabilities at the positions of Super-Kamiokande, 
SNO, and the Gran Sasso Underground Laboratory. 

V. THE ZENITH-ANGLE EXPOSURE FUNCTION AND THE ZENITH-ANGLE 

DISTRIBUTION FUNCTION 

In this section, we define and calculate the zenith-angle exposure function and show how 
the exposure function can be distorted by oscillations and by neutrino-electron interactions 
in the earth into the zenith-angle distribution function. 

The fractional number of neutrino events observed as a function of solar zenith angle, 
a (see Fig. |ID, is, for standard neutrino physics, determined only by the latitude of the 
neutrino detector. The number of events is largest for zenith angles at which the sun 
spends the most time in its apparent motion around the earth. In what follows, we shall 
refer to the normalized number distribution of events, Y{a), as the “zenith-angle exposure 
function”; this function describes the relative amount of time the detector is exposed to the 
sun at a hxed zenith angle. We present in this section the zenith-angle exposure function 
calculated assuming massless neutrinos (no oscillations) and detectors placed at Kamioka, 
Japan (Kamiokande and Super-Kamiokande), Sudbury, Canada (SNO), and the Gran Sasso 
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Underground Laboratory (GALLEX, ICARUS, BOREXINO and HERON/HELLAZ). The 
exposure function for the Homestake detector has been calculated by Cherry and Lande‘ 

m- 


The zenith-angle distribution of events can be distorted by the earth regeneration effect 
implied by MSW solutions, since at solar zenith angles larger than 90*^ neutrinos pass through 
the earth on their way to the detector and therefore z/g’s can be regenerated by interactions 
with electrons in the earth’s interior. The distorted distribution function, /(a), will be 
referred to as the “zenith-angle distribution.” The main goal of this paper is to calculate 
and analyze the shape of /(a) predicted by different MSW solutions and to estimate the 
sensitivity of detectors to the difference between / (a) and the zenith-angle exposure function, 
Y(a). 


A. Zenith-Angle Exposnre Function 

In order to represent accurately the zenith-angle exposure function, we consider 360 
angles, a*, separated by 0.5° intervals between 0 and 180°. The fraction of time in a year 
during which the zenith angle, a, (see Fig. ^ of the sun is close to is proportional to: 

N 

y'{oik) = 9{a{ti) - ak)9{ak+i - a(U)) [1 AU/R(U)]“^ . (4) 

i=l 

Here a{ti) is the solar zenith angle at time ti, N = T/At, T is the duration of one calendar 
year, and At is the time-step. The function 9 is the well known step-function 9{x) = 0, 
X < 0 and 9{x) = 1, a; > 0. The function R{t) is the instantaneous earth-sun distance, and 
1 AU is the one year average earth-sun distance for which solar neutrino fluxes are routinely 
calculated in standard solar models (see Appendix ^ for an explicit representation of the 
time-dependent earth-sun distance). 

The sum in Eq. (^ was computed by simulating the motion of the sun on the celestial 
sphere during one calendar year using the formulae in Appendix Computations with 
different time steps between a few seconds to a few minutes are practically equivalent for 
calculating event rates in the various solar neutrino detectors. 
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The normalized zenith-angle exposure function, Y (ctfc), is obtained directly from Eq. (^, 


N 


r(afc) = F'K)/^F'K). 


(5) 


i=l 


Thus Y[ak) is the fraction of the time during one calendar year that the sun’s zenith angle 
is within the interval {ak — 0.25*^, ak + 0.25°). If the solar neutrino flux is constant in time 
(no neutrino oscillations occur), then the function Y{a) will be the normalized angular 
distribution of events in the detector. We use the analytical expressions for R{t) given in 
Ref. (see Eq. ^2|). The zenith angle exposure function was calculated analytically in 
ref. [^, without including the variation due to the changing earth-sun distance. 

Figures || and |] show the undistorted zenith-angle exposure functions for Super- 
Kamiokande, SNO, and, the detectors assumed to operate at Gran Sasso (ICARUS, BOREX- 
INO, and HERON/HELLAZ). Convenient numerical tables for the Y{ak) are available at 
[http://www.sns.ias.edu/~jnbl . Table |T| lists the latitudes (all northern) of each of the solar 
neutrino detectors. 

The positions of the sharp peaks in Figs. ^ and ^ are determined by the location of the 
detector and by the obliquity, e, of the earth’s orbit (approximately 23.4°). The absolute 
value of the difference between the maximum (or minimum) possible zenith angle for a given 
location and the position of the closest peak is equal to 0 -|- e. At summer solstice, the 
sun’s zenith angle changes from the minimum possible to the angle corresponding to the 
second peak (at an angle > 90°). At winter solstice the sun goes between the peak at angle 
< 90° and the maximum possible zenith angle. Thus during winter solar neutrinos pass 
closer to the earth’s center, whereas during summer they go through lower density layers of 
the mantle. 


B. Zenith-Angle Distribution Function 

We calculate the event rates, Qi, along each direction, with the aid of a set of survival 
probabilities, Pse: computed (just once for each direction, mixing angle, and E/Am'^) along 
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a fixed set of trajectories through the earth. The numerical value of each Qi is obtained 
from Eq. E3 of Appendix C. The normalized zenith-angle distribution function is 

f{ai) = Q{ai)Y{ai)/^Q{ak)Y{ak). ( 6 ) 

k 

In the absence of oscillations, the detector event rate, Qi, is independent of direction and 
disappears from the right-hand side of the equation. In this case /(a) = Y{a). 

In Figs. I and H, we present the expected distorted angular distributions for the 
SMA, LMA, and LOW solutions in Super-Kamiokande, SNO, ICARUS, BOREXINO and 
HERON/HELLAZ. In the panels for Super-Kamiokande, SNO and ICARUS, we show only 
the distributions corresponding to the best-ht points in the SMA and LMA solutions; the 
curves corresponding to the LOW solution are virtually undistorted. Correspondingly, in 
the panels for BOREXINO and HERON/HELLAZ, we show only the expected angular dis¬ 
tribution for the LOW solution (see Eq. since the SMA and LMA solutions imply only 
a negligible distortion of the zenith-angle distribution for these low energy detectors. 


VI. MOMENTS OF THE ZENITH-ANGLE DISTRIBUTION 


In this section, we evaluate the expected MSW distortions of the angular distribution 
and compare the hrst two moments of the predicted angular distribution of events with the 
calculated moments expected in the absence of oscillations. This comparison constitutes a 
new and, for the small mixing angle solution, a more powerful way of analyzing the time 
dependence of the observed neutrino events. The predicted distortions of the recoil electron 


energy spectra in Super-Kamiokande and in SNO were investigated in Ref. |]^,|3^ in terms 
of the analogous moments of the energy distribution. 

In Sec. |V1 A| , we dehne the hrst and second moments of the zenith-angle distribution and 
then calculate in Sec. |VlB| the predicted MSW changes in the zenith-angle distributions. We 
plot the relative shift of the hrst moment in the plane of the neutrino oscillation parameters 
Am^ and sin^ 29, illustrating the range of possible values of the shift of the hrst moment 
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within the 95% C.L. allowed by the Chlorine |l|, Kamiokande [0, GALLEX ||^, and SAGE 
^ experiments. 

How sensitive will the Super-Kamiokande, SNO, IGARUS, BOREXINO and HEL- 
LAZ/HERON experiments be to the regeneration effect? We answer this question in 
Sec. |V1 B| by plotting the number of standard deviations each MSW solution is separated 
from the no-oscillation solution in the plane defined by the values of the first two moments. 


A. Moment Definitions 


The first two moments of the zenith-angle distribution are defined by 


(a) = j af{a)da, 


(7a) 


and 


= y~ {a)Yf{a)da. 


(7b) 


The fractional shifts in these two moments are defined as 


and 


tto 


((a) — ao)/ ao, 


(8a) 





(8b) 


where Oq and cTq are the moments of the zenith-angle exposure function distribution (no 
oscillations) and (a) and are the moments of the distorted (by oscillations) zenith-angle 
distribution. Because of the symmetry of Y{a) about a = 90 deg, ao = n/2 for all the 
detector locations. The values of the undistorted second moment, cxo, are: 37.9° (Super- 
Kamiokande), 32.9° (SNO), 34.9° (Gran Sasso), 33.9° (Homestake), 34.5° (Baksan), and 
47.5° (Equator). 

The Kamiokande collaboration excluded a significant region in Arri^ - sin^ 29 parameter 
space by grouping the events from below the horizon into five bins, corresponding to different 
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zenith angles P]. Lisi and Montanino in Ref. 1^^ calculated binned angular distributions 
for Super-Kamiokande and SNO for different values of Am^ and 9“^. 


The calculation of the hrst two moments of the zenith-angle distribution is subject to 
fewer complications than the more traditional method used in references P,|3^ of binning the 
data and doing a analysis. For binned data, the a priori unknown normalization and the 
angular dependence of the detector sensitivity are, for example, two of the aspects of neutrino 
experiments that introduce bin-to-bin correlations which are often difficult to estimate from 
the observations or to evaluate including all of the relevant detector characteristics. On 
the other hand, the calculation of the first and second moments directly from the data is 
st raight forward. 


B. Predicted MSW Moments 

Figures ^ and ^ show the expected shifts in the first two moments as a function of the 
neutrino oscillation parameters. Am? and sin^ 29. We plot in Fig. (Fig. ^ contours of con¬ 
stant fractional shift, Aa/ao (Acr^/cTg), of the first moment (second moment) in the Am^ - 
sin^ 20 plane. The five panels are, for both figures, for Super-Kamiokande, SNO, ICARUS, 
BOREXINO and HERON/HELLAZ. For Super-Kamiokande, Aa/ag varies between —0.1% 
and 3% in the SMA region and between 0.5% and 7% in the LMA region. The corresponding 
ranges for SNO are (—0.1,4.5)% (SMA) and (0.3,15)% (LMA). SNO is somewhat more sen¬ 
sitive than Super-Kamiokande to the shift of the first moment because z/^s (and Urs) do not 
contribute to the charged-current signal in SNO. ICARUS and SNO have similar sensitivi¬ 
ties to the regeneration effect since both observe charged current reactions with high-energy 
thresholds of 10.9 MeV and 6 MeV, respectively. The low-energy experiments, BOREXINO 
and HERON/HELLAZ, are insensitive to both the SMA and LMA allowed regions, but will 
be able to test the LOW solution to which neither SNO nor Super-Kamiokande are sensitive. 
The results for the second moment, shown in Fig. exhibit similar trends. 

Figure ^summarizes the potential of the second generation of solar neutrino experiments 
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for discovering new physics via the earth regeneration effect. The hgure displays iso-sigma 
ellipses, statistical errors only, in the plane of the fractional percentage shifts of the hrst 
two moments, Aa/cto and Acr^/cTg. Assuming a total number of events of 30000, (which 
corresponds to ~ 5 years of standard operation for Super-Kamiokande and ~ 10 years 
for SNO), we have computed the sampling errors on the hrst two moments as well as the 


correlation of the errors using the following well known formulae cr(mi) = A = 


, and p(mi, m 2 ) =-Here Ui is the i—th moment of the angular distribution 

\/M2VM4-M2 

of events and m* are the corresponding estimates of the moments from a sample of size 
N (number of events). The iso-sigma ellipses for the six detectors we consider here are 
centered around the undistorted zenith-angle exposure function for which, by dehnition, 
Aa = A(T^ = 0. Figure ^ shows for each detector the predicted shifts of the hrst two 
moments in the SMA, LMA, and LOW solutions. The horizontal and vertical error-bars 
denote the spread in predicted values of the shifts in the hrst two moments, which are 
obtained by varying Am^ and sin^ 29 within the 95% C.L. allowed (see Fig. |^) by the four 
pioneering solar neutrino experiments. 

For Super-Kamiokande (SNO), the current best-ht parameters Am^ and sin^26' [see 
Eq. (j^)] predict a 5 ct (6.5ct) ehect for the SMA solution and 13 ct (25(t) ehect for the LMA 
solution. Note that SNO is expected to require twice as much time to collect the same 
number of events as Super-Kamiokande. In the same amount of observing time, SNO and 
Super-Kamiokande are approximately equivalent for the SMA and Super-Kamiokande is 
signihcantly more efficient for the LMA. 


VII. DAY-NIGHT ASYMMETRY 


In this section, we calculate the day-night asymmetry |1l3|- [T5| , |3T|j33| -p5| caused by the 
earth regeneration ehect. For Super-Kamiokande and SNO, the hrst accurate calculations of 
the day-night asymmetry have been made only recently using a realistic density distribution 
in the earth 
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We evaluate numerically the day-night asymmetry for Super-Kamiokande, SNO, 
ICARUS, BOREXINO, and HERON/HELLAZ. We present quantitative estimates of the 
sensitivity of future experiments to the predicted difference in night-time and day-time event 
rates. 

The day-night asymmetry is dehned as: 




Qn Qd 


(9) 


Qn + Qd^ 

where Qn and Qd are respectively the average night-time and day-time event rates during 
one year. The calculation of the asymmetry is made using the one-year average night-time 
and day-time survival probabilities described in Appendix A. 

Figure shows contours of constant day-night asymmetry in the Am^ - sin^ 26 plane 
for Super-Kamiokande, SNO, ICARUS, BOREXINO, and HERON/HELLAZ. The shaded 
regions in the panels for Super-Kamiokande, SNO, and ICARUS are allowed (at 95% C.L.) by 
the present data from Homestake, Kamiokande, GALLEX and SAGE. In order to include 
the LOW solution near Am^ = 10“^ eV^, the panels corresponding to BOREXINO and 
HERON/HELLAZ show the allowed regions at 99% C.L. 

Table gives the expected range of day-night asymmetries in Super-Kamiokande and 
SNO for both the SMA and LMA solutions. The ICARUS sensitivity is similar to SNO. 
In general, SNO is slightly more sensitive than Super-Kamiokande for testing the day-night 
asymmetry because SNO detects a CC reaction; the sensitivity to NC scattering of muon 
neutrinos decreases somewhat the sensitivity of the Super-Kamiokande detector. Both de¬ 
tectors are very sensitive to the large day-night asymmetry predicted in the LMA solu¬ 
tion. Their sensitivity to the SMA solution is smaller, especially in the low sin^ 29 corner 
(sin^20 < 5 X 10“^), where An-d — 0. Both Super-Kamiokande and SNO are practically 
insensitive to the LOW solution because they can register only high-energy boron neutrinos 
which do not resonate in the earth for the low Am^ values in the LOW solution. As noted in 
Ref. [M, An-d is negative in a small region of parameter space for both Super-Kamiokande 
and SNO, i.e., the night time event rate can be slightly lower than the day-time event rate. 
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The differences between SNO and ICARUS are mainly due to the different assumed neu¬ 
trino thresholds (6.4 MeV and 10.9 MeV, respectively) and to the locations of the detectors. 
For large mixing angles, for which the regeneration effect takes place mainly in the man¬ 
tle, the sensitivities of the two detectors are particularly similar. At small mixing angles 
(sin^ 26* < 0.3) and 5 x 10“® < Am^/eV^ < 10“^, ICARUS is slightly more sensitive to the 
earth effect because it is located at a lower latitude than SNO. 

The three future detectors that will measure low-energy neutrinos, BOREXINO (^Be), 
HERON {pp), and HELLAZ {pp), will be sensitive to the large day-night asymmetries pre¬ 
dicted in the LOW solution, but insensitive to the asymmetries predicted by the SMA and 
LMA solutions. The range of day-night asymmetries expected for the LOW solution are 
(16 ± 5%) for HERON/HELLAZ and (16 ± 8%) for BOREXINO. 

VIII. WHICH STATISTICAL TEST IS BEST? 

Which statistical tests are most powerful in detecting new physics? What type of analysis 
will most clearly show departures from the zenith-angle exposure function due to the regen¬ 
eration effect? By analyzing simulated data in this section, we shall see that the preferred 
statistical analysis depends upon which solution Nature has chosen. 

Table 0 compares the sensitivity of Super-Kamiokande and SNO to the earth regener¬ 
ation effect for three different statistical tests. We have computed the number of standard 
deviations by which the best-ht MSW solutions (described in Sec. ^V]) differ from the undis¬ 
torted zenith-angle exposure function. We consider the hrst and second moments of the 
zenith-angle distribution (see Section 0), the day-night asymmetry {A^-d) (see Sec. 0). 
and the Kolmogorov-Smirnov test of the distorted zenith-angle distribution. We assume 
30000 events are detected in the case of the SMA solution. The comparison is made after 
only 5000 events are observed for the more-easily recognized LMA solution. 

For the SMA solution, the moments analysis is most sensitive. The difference for SNO, 
between 4.9a (day-night analysis) and 6.5a (moments distribution), corresponds to 20000 
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events, or approximately 7 years of data taking. All three statistical tests can easily reveal 
the best-ht LMA solution, although the day-night asymmetry is the most efficient char¬ 
acterization in this case. The Kolmogorov-Smirnov test is the least sensitive to the SMA 
solution, but performs better than the moments method for the LMA solution. 

We can understand physically why the SMA distortion is most easily detected by measur¬ 
ing the moments while the LMA distortion is most prominent in the day-night asymmetry. 
Figure 0 shows for Super-Kamiokande the fractional distortion, [/(a) — Y{a)]/Y{a), of 
the zenith-angle distribution for the best-ht SMA and LMA solutions. One can crudely 
approximate the distortions by, for the SMA solution, a delta-function near the maximum 
allowed zenith angle and, for the LMA solution, a step function near 7r/2. With these 
simple approximations, one can show analytically that the hrst moment and the day-night 
asymmetry have similar statistical power for the SMA solution, and the second moment is 
more discriminatory than either the hrst moment or the day-night asymmetry. The reason 
that the second moment is so useful for the SMA solution is that in this case the distortion 
mostly arises when the neutrinos pass through the core at large zenith angles. Because the 
vacuum mixing angle is small, the enhanced mixing [p^es = 7 g cm“^(i7/10 MeV)] due to 
the earth matter ehect is particularly signihcant when the neutrinos traverse the core. Since 
the vacuum mixing angle is large for the LMA, the matter enhancement is not especially 
signihcant in this case and the main regeneration for the LMA solution is due to oscillations 
that occur in the mantle, i.e., whenever a > 7r/2. The day-night asymmetry is well-tuned 
to this distortion since An_d compares the average event rate for a > 7r/2 with the event 
rate for a < tt/2. The LOW solution produces a relative distortion, [/(a) — Y{a)\/Y{a), 
that has a shape similar to the LMA solution and is therefore most easily detected by the 
day-night asymmetry. The Kolmogorov-Smirnov test is not optimally tuned to any of the 
three best MSW solutions and is therefore not as powerful as the moments or the day-night 
asymmetry. 
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IX. MOMENTS OF THE ENERGY SPECTRUM 


We refine in this section our previous calculations of the hrst two moments of the energy 
spectrum from electron recoils produced by interactions with neutrinos. We use here 
the slightly improved MSW solutions, described in Sec. that include regeneration in 
the earth. The reader is referred to our earlier paper for the relevant dehnitions and 
notation (see also ref. ||^] for a similar calculation). 

Table |V11| presents, for an assumed threshold of 5 MeV, the hrst and second moments of 
the electron recoil energy spectrum, and the percentage shifts with respect to the average 
electron kinetic energy, Tq, and the dispersion in the kinetic energy, do, in the absence of 
oscillations. The results for Tq and ao differ by small amounts (< 1%) from our earlier 


results given in ref. [50|; the present results are numerically more precise. The results in 


Table |VII| are given for the best-estimate MSW solutions (SMA, LMA, and LOW) described 
in Sec. B For completeness, we list the one year average moments of the energy spectrum for 
day-time, night-time, and the total year. The calculated moments for the Super-Kamiokande 
experiment are given in the upper part of the table and the moments for SNO are listed in 
the lower part of the table. Table |VIII| presents the same results for an assumed threshold 
of 6 MeV. 

We did not include in our calculation the unknown trigger efficiencies of Super- 
Kamiokande or SNO. The inclusion of these trigger functions can change the predicted 
hrst and second moments of the energy distribution by a few percent and will certainly be 
included in the careful Monte Carlo calculations that will be performed ultimately by the 
Super-Kamiokande and SNO experimental groups. 


Comparing the calculated day and night rates. Table |VI1| and Table |VIII| show that 
regeneration in the earth slightly decreases , for both the SMA and the LOW solutions, 
the average kinetic energy of the recoil electrons in both Super-Kamiokande and SNO. This 
decrease occurs because in the sun these two solutions preferentially transform low energy 
neutrinos from i/g to (or u-r) and therefore there is a relatively larger chance at low energy 
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of regenerating z/g from z/^ ( or i>r) in the earth. For the LMA solution, regeneration increases 
the average kinetic energy since in this case the high-energy part of the neutrino energy 
spectrum is preferentially depleted of z/g in the sun. 

The shift between day and night of the moments is most signihcant for the LMA solution. 
In fact, if the nature has chosen the LMA solution, then the spectral distortion may be 
highlighted by comparing the day-time and night-time moments. 


X. SENSITIVITY TO EARTH MODELS AND SOLAR MODELS 


We calculate in Sec. pC A| the sensitivity of the MSW predictions to the assumed density 
prohle and chemical composition of the earth model and in Sec. PCB| the dependence upon 
the assumed model of the sun. 


A. Uncertainties Due to Earth Models 


Table | presents the calculated percentage shifts of the hrst two moments of the zenith- 
angle event distribution for all six models of the earth discussed in Sec. B the calculations 
were made assuming either the SMA or the LMA solutions. The fractional changes of the 
hrst moment vary by only ~ 0.02% for the SMA solution and ~ 0.2% for the LMA solution, 
although the density prohles in some of these models are signihcantly different from the 
range allowed by current seismological data. We conclude that the shape of the zenith-angle 
distribution can be calculated with acceptable accuracy for any of the recently published 
density prohles of the earth. 

Table illustrates the uncertainties in the MSW predictions of the shifts of the hrst 
and second moments of the zenith-angle distribution due to uncertainties in the electron 
number density in the mantle and in the core. The ranges of Z/A included in the table 
(±2% in the core and —1%, —2% in the mantle) are larger than the current estimates of the 
geophysical uncertainties (see the discussion in Sec. B and reference P^,^). We conclude 
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from Table |I^ that uncertainties in the chemical composition affect the predicted moment 
shifts due to regeneration by at most a few percent of their values. 

A simplihed model with a uniform composition of ZjA = 0.5 has been used in refer¬ 
ence (and in many or most of the early calculations related to the regeneration effect, 
see 1^). The predictions from this constant-composition model are also given in Table 
this crude model leads to imprecise, but not grossly erroneous, predictions of the moments 
of the zenith-angle distribution. 


B. Uncertainties Dne to Solar Models 


To the best of our knowledge, all previous discussions of the earth regeneration effect 
have described this phenomenon as if it were completely independent of solar models. This 
implicit assumption is not exactly correct since the size of the earth regeneration effect 
depends upon the flavor content of the incident neutrino beam, which must be calculated 
by using a solar model to describe (for specihed MSW parameters) the production and 
conversion probabilities of solar neutrinos as a function of the position in the sun at 
which the neutrinos are created and the neutrino energy. The slightly different density 
distributions in different solar models have the largest effect, which is still quite small as we 
shall see below, on the inferred flavor content of the incident solar neutrino flux. 

In order to quantify the dependence of the predicted earth regeneration effect upon the 
characteristics of the solar model, we have calculated the fractional shifts of the hrst and 
second moments of the zenith-angle event distribution using three different solar models. 
As our standard solar model, we adopt the model with helium and heavy element diffusion 
of Bahcall and Pinsonneault P]. For comparison, we use the 1992 model of Bahcall and 
Pinsonneault [^, which includes helium diffusion (but not heavy element diffusion) and 
somewhat less accurate input physics. Finally, we use the 1988 model of Bahcall and Ul¬ 
rich |5^, which does not include any diffusion and has less precise opacities, equation of 
state, and other input data. 
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Table shows that the MSW predictions are essentially identical for the 1992 solar 


model with helium diffusion and the 1995 solar model with helium and heavy element diffu¬ 
sion plus improved input data. The 1988 solar model leads to predictions that can differ by 
as much as 10% for the SMA moments that will be measured by SNO. However, this 1988 
model is inconsistent with recent helioseismological measurements since the 1988 model does 
not include diffusion [^. 

We conclude that predictions of the earth regeneration effect are practically independent 
of solar models as long as the models include diffusion (i.e., are consistent with helioseismol¬ 
ogy) • 


XI. FUTURE EXPERIMENTS AT THE EQUATOR 

Recently Gelb, Kwong and Rosen suggested building a new detector similar to SNO 
close to the equator in order to increase the sensitivity of the experiment to the earth 
regeneration effect. An equatorial location maximizes the time neutrinos pass through the 
core of the earth during one calendar year. 

In this section, we calculate the size of the regeneration effect for hypothetical equa¬ 
torial detectors and compare with the sensitivity of the detectors in their actual posi¬ 
tions. We consider equatorial analogues of the Super-Kamiokande, SNO, BOREXINO, and 
HERON/HELLAZ detectors. 

Figure |T^ shows the predicted zenith-angle distribution for detectors at the equator. The 
curve in the upper left panel is the zenith-angle exposure function. The other Eve panels 
show the distortion due to regeneration for equatorial analogues of Super-Kamiokande, SNO, 
ICARUS, BOREXINO, and HERON/HELLAZ. For each detector, only the predicted angu¬ 
lar distribution functions are shown for the best-ht solutions to which the relevant detector is 
sensitive: SMA and LMA for the high-energy boron neutrino detectors (Super-Kamiokande, 
SNO and ICARUS) and LOW for the low-energy neutrino detectors (BOREXINO and 
HERON/HELLAZ). 
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Figure |T^ shows the iso-sigma ellipses for the four equatorial detectors. BOREXINO 
and HERON/HELLAX would be more sensitive to the LOW solution if these detectors 
were built at the equator. However, the high-energy neutrino detectors (Super-Kamiokande, 
SNO, and ICARUS) would remain insensitive to the LOW solution even if they were moved 
to the equator. Low-energy neutrino detectors (BOREXINO and HERON/HELLAZ) also 
remain insensitive to the SMA and LMA solutions even at the equator. 

Table ^ shows the gain in sensitivity that would occur if detectors like Super-Kamiokande 
and SNO were built at the equator. The enhancement is represented in the table by the 
shift in the hrst and second moment of the zenith-angle distribution and by the day-night 
asymmetry. The enhancements would be important for the best-ht SMA solution, but less 
signihcant for the LMA solution. However, regions of parameter space in the LMA solution 
for which the predicted shifts in the hrst and second moments are small could be probed 
more precisely with detectors at the equator. 

XII. DISCUSSION AND CONCLUSIONS 

The conversion in the earth of (or to the more easily detected Uf. is a distinctive 
prediction of the MSW effect that offers the possibility of unambiguously establishing the 
existence of physics beyond the standard electroweak model. Because of the importance of 
this subject, we have carried out precise numerical calculations of the size of the regeneration 
effect predicted by different MSW parameters that are consistent with the experimental 
results from the chlorine, Kamiokande, GALLEX, and SAGE experiments. Our results show 
the potential of the new experiments, Super-Kamiokande, SNO, IGARUS, BOREXINO, 
HERON, and HELLAZ, for discovering the regeneration effect. 

Our results provide the most precise predictions available of the expected zenith-angle 
distribution of the solar neutrino events in the absence of new physics and in the presence of 
MSW distortions. The results are obtained by numerical calculations that are discussed in 
Sec. |V| and illustrated in Fig. ^ (for Super-Kamiokande and SNO) and Fig. ^ (for the Gran 
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Sasso experiments ICARUS, BOREXINO, and HERON/HELLAZ). 

We present the predictions for the small mixing angle (SMA), large mixing angle(LMA), 
and low mass (low probability, LOW) MSW solutions of the solar neutrino problems. The 
parameters of these MSW solutions, which are consistent with the results of the chlorine, 
Kamiokande, GALLEX, and SAGE experiments, are given in Sec. Our solutions include 
self-consistently the effects of earth regeneration. 

Figure § shows the allowed regions of the three MSW solutions in the Am^ - sin^ 29 plane. 
Figure | presents the survival probabilities as a function of energy for created in the sun. 
This figure compares survival probabilities computed for the day (without regeneration) 
with survival probabilities for the night (with regeneration) and with the average annual 
survival probabilities. 

We describe the predicted MSW distortions in terms of the first two moments of the 


zenith angle distribution of neutrino events (see Sec. [V^ , as well as in terms of the traditional 
day-night asymmetry (see Sec. ivg. We analyze simulated data in Sec. |V111| and show that 
the moments of the zenith-angle distribution are more sensitive to the harder-to-detect SMA 
solution. The predicted large effect of the LMA solution is more easily discovered with the 
conventional day-night asymmetry. 

The “bottom line” is illustrated succinctly in Fig. This figure shows that the current 
best-estimate MSW solutions predict statistically signihcant deviations from the undistorted 
zenith-angle moments for the Super-Kamiokande, SNO, and IGARUS experiments (which 
are sensitive to the SMA and LMA solutions) and the BOREXINO and HERON/HELLAZ 
experiments (which are sensitive to the LOW solution). 

We have considered a number of effects that have not been previously investigated in 
connection with the earth regeneration effect. We have calculated the sensitivity of the 
MSW predictions to a wide range of density prohles of the earth and also to a set of extreme 
chemical compositions. These calculations are discussed in Sec. 0 and Sec. pC A[ We also 
evaluate the slight dependence of the predicted earth regeneration effect upon the assumed 
solar model used to calculate the flavor content of the incident neutrino beam (see Sec. |XB| ). 
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Our results show that these usually-neglected effects associated with the earth and solar 
models are rather small. 


For completeness, we have carried out calculations for hypothetical new detectors that 
might be built near the equator. These calculations are described in Sec. 0 and show quan¬ 
titatively the enhanced sensitivity to the earth regeneration effect of equatorial detectors, 


as emphasized by Gelb et al. |]28 


Using the best-fit MSW solutions calculated here that include the earth regeneration 
effect, we have evaluated the first and second moments of the electron recoil energy spectrum 
for neutrinos detected in Super-Kamiokande and SNO. These calculations, summarized 
in Sec. ^ and in Table |VII| and Table [VII1| , refine our earlier results (see also ref. p4| ) 
for the moments of the electron recoil spectrum. Perhaps most importantly, they show that 
for the LMA solution the comparison of the recoil electron energy spectrum between day 
and night may reveal a distortion that is not apparent in the temporal average of the energy 
spectrum. 

What would we learn from an observation which showed that the neutrino counting rate 
depended upon solar direction? The experimental demonstration of a dependence of solar 
neutrino event rate upon the direction of the sun would not only constitute a direct proof of 
new physics but would at the same time eliminate a number of the popular alternatives to 
the MSW effect. Many of the alternatives to the MSW effect, such as vacuum oscillations, 
magnetic moment transitions, and violations of the equivalence principle predict that the 
counting rate is independent of the zenith angle position of the sun. 
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APPENDIX A: NEUTRINO SURVIVAL PROBABILITIES 


In order to calculate the event rates as a function of time we hrst compute, following the 
prescription in Ref. [^, the electron neutrino survival probabilities, Pse, after traversing 
the earth. We begin by using the analytical approximation developed in Ref. of the 
survival probabilities, Ps, for an electron neutrino passing through the sun; these solar 
survival probabilities are averaged over the relevant neutrino production regions for ®B, ^Be, 
pp, pep, and CNO neutrinos. Assuming that the neutrinos arriving at the earth represent an 
incoherent superposition of mass-eigenstates we calculate the electron neutrino survival 
probability after passing through the earth, Pse, from the expression: 


Pse — 


Ps - sin^ 9 + P 2 e(l - 2Pg 
cos 26 


(Al) 


Here 6 is the neutrino mixing angle in vacuum and Ps is the probability of an electron 
neutrino produced in the sun to arrive as an electron neutrino at the earth. The ambiguity in 
the above equation for maximal mixing (cos 26 = 0) is only apparent. It can be shown that in 
this case both Ps and Pse are equal to 1/2. The probability P 2 e to hnd an electron neutrino 
after passing through the earth, if the initial state describing a neutrino entering the earth is 
the pure mass-eigenstate 1 ^ 2 , is calculated numerically by integrating the evolution equations 
for neutrino states in matter in the form given in A key point in the calculation is that 
these survival probabilities as a function of E//S.w? need be calculated just once for each 
sin^ 26 and for a hxed set of trajectories. We have chosen a set of 180 trajectories equally 
spaced between 90 and 180 degrees of zenith angle separated by 0.5 degrees. This set is 
the same for each solar neutrino detector since the density prohles used in our calculations 
and described in Sec. El are spherically symmetric. The survival probabilities along each 
trajectory, Pse, have been calculated for a grid of 201x201 values of E//S.w? and sin^ 26\ 


10^ MeV/eV^ < P/Am^ < 10^^ MeV/eV^ 


(A2a) 
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and 


10“^ < sin^ 29 < 1. (A2b) 

The numerical precision of the calculated survival probabilities along each trajectory is better 
than 0.1%. 

The one-year averaged survival probability is given by: 

N 

PsE^T.P'sEY{ai), (A3) 

i=l 

where the sum is over zenith angles from 0 to 180 degrees. (iV = 360 in our case). The 
zenith-angle exposure function is defined in Sec. Correspondingly, the one-year averaged 
night-time survival probability is given by 

N 

PSe= E P!iEy(ai), (A4) 

i=Nf2 

where the sum now runs over angles from 90 to 180 degrees. Since the night-time and 
day-time intervals within one year are equal, the day-time event rate is simply 

= O.bPsE. (AS) 

With the calculated survival probabilities, it is straightforward to determine the correspond¬ 
ing total event rates, as well as the day-time and night-time one-year averaged event rates 
in any solar neutrino detector. 


APPENDIX B: TIME DEPENDENCE OF THE ZENITH ANGLE 


The dependence of the solar zenith angle on the time of the year and on the geographic 
location of the detectors is given by the following set of formulae ||49|| : 

cos a = sin 5 sin 0 -|- cos 6 cos cj) cos H , (Bla) 


sin 6 = sin e sin A , 


L = 280°.461 + 0°.9856003n , 


(Bib) 

(Blc) 


29 




n = -1462.5+ L> + i/, 


(Bid) 


g = 357°.528 + 0°.9856003n, (Ble) 

and 

\ = L + 1°.915 sin g + 0°.020 sin 2g . (Blf) 

The precision in the apparent coordinates of the Snn is 0°.01 and the precision of the 
eqnation of time is 6 seconds between the years 1950 and 2050. Here H is the fraction of 
day from 0^ UT, D is day of the year (connting from Jannary 1), n is the nnmber of days 
from Jnlian year 2000.0, A is the ecliptic longitnde, L is the mean longitnde of the snn 
(corrected for aberration), e = 23°.439 — 0.0000004n is the obliqnity of the ecliptic, 6 is the 
snn’s declination, and g is the mean anomaly]^. 

The distance of the snn from the earth in astronomical nnits (1 AU = 1.495978706(2) 10^^ 
m) is given by the formnla: 

R = 1.0014 - 0.01671 cos g - 0.00014 cos 2g. (B2) 

Eqnation (p^) has been nsed in the calcnlation of the day-night asymmetries and the shifts 
of the first two moments of the zenith-angle distribntion. 


APPENDIX C: CHARACTERISTICS OF FUTURE DETECTORS 


We describe in this Appendix the characteristics we have assnmed for Snper-Kamiokande 
SNO O, ICARUS [0], BOREXINO |^, HERON and HELLAZ |^. The Snper- 


Kamiokande, SNO and ICARUS detectors are sensitive only to high-energy ®B nentrinos, 
while BOREXINO is sensitive primarily to ^Be nentrinos {E^ = 0.862 MeV) and the HERON 
and HELLAZ detectors are being designed to detect low energy {Ey < 0.44 MeV) pp nen¬ 
trinos. 


For definitions of these astronomical qnantities see Ref. 
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For the Super-Kamiokande detector, we adopt a threshold of 5 MeV and a trigger ef- 


hciency of 50% at this energy [^. The energy resolution function is assumed to have a 
gaussian shape with FWHM of 1.6 MeV at electron kinetic energy 10 MeV. We have per¬ 
formed calculations, which show that the sensitivity of our results to the assumed energy 
resolution and trigger efficiency of the detector. 

For the SNO detector, we calculate only the rate of the CC reaction, namely + d ^ 
p + p + e~. We adopt a threshold of 5 MeV and an energy resolution function with a la 

uncertainty of 1.1 MeV at 10 MeV electron kinetic energy. The CC cross-section for SNO 


was taken from Ref. . The trigger efficiency function has been approximated with a step 
function at the threshold of the detector. 

Reference gives further details regarding our characterization of SNO and Super- 
Kamiokande. 

For ICARUS, we have considered only the superallowed transition and have used the 
neutrino absorption cross sections given in ref. [Q. We have assumed a neutrino threshold 
for detection that corresponds to electrons being produced with at least 5 MeV of energy, 
which requires a minimum neutrino energy of 10.9 MeV. 

The BOREXINO detector is being developed as a neutrino-electron scattering exper¬ 
iment that will measure the flux of ^Be neutrinos, using the step in the energy distribu¬ 
tion of neutrino-electron scattering events at the maximum recoil-electron kinetic energy of 
Te = 0.62 MeV. The detector characteristics that are knowable a priori (which does not 
include the crucial background rate versus energy) are not as important as for SNO and 
Super-Kamiokande and we need to compute only the z/g survival probability. We include 
radiative corrections to the neutrino-electron cross-section, calculated in |^, which for re¬ 
coil electron energies below 0.62 MeV are less than 1%. The ratio of the {r'^e) to (z/gc) total 
cross-sections for neutrino energy 0.862 MeV is a^i^^e/o'ue,e = 0.221. 

HERON and HELLAZ are also neutrino-electron scattering experiments with very dif¬ 
ferent preliminary designs, but with the same target material, helium. They will measure 
the flux and spectral shape of pp neutrinos. We have assumed a threshold of 0.1 MeV 
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and a perfect energy resolntion (delta fnnction). The trigger efficiency is represented by a 
step-fnnction at the threshold of the detector. We again nse nentrino-electron cross-sections 
inclnding radiative corrections ^ . 


The event rate, Q, averaged over certain time interval, r, in a nentrino-electron-scattering 
experiment, snch as Snper-Kamiokande, BOREXINO, HERON or HELLAZ is given by 


{Q)r = I [Z,SE.){Pm)r + (1 - {Pm)r) 


dE„ 


(Cl) 


Here are the response fnnctions of the detector to either Ve or z/^, <h is the solar 

nentrino flnx to which the detector is sensitive, Ey is the nentrino energy, E^ax is the 
end point of the nentrino energy spectrnm and {P)t is the average snrvival probability 
for the chosen time interval r. The interval can be, e.g., the total day-time or night-time 
dnring one calendar year, or a whole calendar year inclnding both days and nights. The 
calcnlation of the average snrvival probabilities is described in Appendix ^ (Eqs. |A1[ - |A5D . 
The response fnnctions represent the convolntion of the absorption cross sections with the 
detector characteristics (see Ref. |OT for details). 
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TABLE 1. Sensitivity to the model of the earth. The table illustrates the weak dependence on 
the model of the earth of the calculated changes in the first and second moments of the angular 
distribution of events in the Super-Kamiokande and SNO detectors. The density distributions in 
the six models of the earth listed in the table span a range of possibilities that is much larger than 
suggested by current geophysical knowledge. The second and third columns give the total mass, 
M 0 (in 10^'^ g), and the moment of inertia, I (in lO^^gcm^), for each model. The mass of the earth 
is 5.97370 ± 0.00076, the polar value of the moment of inertia is 0.804, and the equatorial value 
is 0.801 1^.] The last four columns give the fractional shift in percent of the first two moments 
(Aa/ao and A(T^/(Tq) of the zenith-angle distribution of events in the Super-Kamiokande and SNO 
detectors for the SMA (upper entry) and for the LMA (lower entry) solutions. 


A{a)/ao 

Model M 0 I Super-Kamiokande (%) 

AuVcrg 

Super-Kamiokande (%) 

A{a)/ao 

SNO (%) 

Ao-Vf^o 

SNO (%) 

PREM 5.972 0.802 1.02 

2.03 

1.04 

2.42 
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1981 


3.23 


-0.33 


5.28 


- 2.20 


hb; 

5.970 

0.801 

1.00 

2.00 

1.02 

2.38 

1975 



3.29 

-0.27 

5.37 

-2.09 

B497 

5.067 

0.800 

1.01 

2.03 

1.03 

2.41 

1973 



3.22 

-0.26 

5.26 

-2.05 

B1 

5.949 

0.798 

1.01 

2.03 

1.03 

2.42 

1974 



3.23 

-0.33 

5.21 

-2.38 

A" 

5.946 

0.798 

1.00 

2.01 

1.01 

2.39 

1967 



3.29 

-0.51 

5.35 

-2.55 

Stacey 

5.973 

0.802 

1.00 

2.01 

1.02 

2.38 

1969 



3.27 

-0.28 

5.34 

-2.11 


TABLE II. Solar neutrino data used in the analysis. The experimental results are given in SNU 
for all of the experiments except Kamiokande, for which the result is expressed as the measured 
flux above 7.5 MeV in units of cm“^s“^ at the earth. The ratios of the measured values to the 
corresponding predictions in the standard solar model of Ref. are also given. The result cited 
for the Kamiokande experiment assumes that the shape of the ®B neutrino spectrum is not affected 
by physics beyond the standard electroweak model. 


Experiment 

Result 

Theory 

Units 

Result/Theory 

Reference 

HOMESTAKE 

2.56 ±0.16 ±0.14 


SNU 

0.27 ±0.022 

i 

KAMIOKANDE 

2.80 ±0.19 ±0.35 

6.62 

10 ® cm“^ s“^ 

0.42 ± 0.060 

0 

GALLEX 

69.7 ±6.7 ± 

136.81^ 

SNU 

0.51 ±0.058 

i 

SAGE 

rjcy +12 +5 
-10 -7 

136.81^ 

SNU 

0.53 ±0.095 

i 
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TABLE III. Northern latitudes (in degrees) for several solar neutrino detectors. These angles 


correspond to 0 in Fig. ||. 


Homestake 

Kamioka 

Gran Sasso 

Baksan 

Sudbury 

44.3 

36.4 

42.4 

43.3 

46.3 


TABLE IV. Sensitivity to Solar Models. Values of the percentage shifts of the first and second 
moments of the zenith-angle event distribution in Super-Kamiokande and SNO computed with 
three different solar models (BU1988, BP1992, BP1995). 




SMA 


LMA 


Detector 

Solar Model 

A(a)/ao 

AcrVc^O 

A(q;)/q;o 

AdVc^o 


1995 

1.02 

2.03 

3.22 

-0.32 

Super-Kamiokande 

1992 

1.02 

2.03 

3.22 

-0.32 


1988 

0.96 

1.94 

3.22 

-0.32 


1995 

1.04 

2.42 

5.28 

-2.20 

SNO 

1992 

1.04 

2.44 

5.28 

-2.20 


1988 

0.96 

2.26 

5.27 

-2.20 
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TABLE V. Day-night asymmetry in Super-Kamiokande and SNO. The table gives the magnitude 


of the expected day-night asymmetry [A^-d, see Eq. (|9|)) (in percent) in Super-Kamiokande and 
SNO for values of the neutrino oscillation parameters Am^ and sin^ 29 corresponding to the best-fit 


SMA and LMA solutions (see Eqs. ]^,1T and 2a ,pb|) . The indicated uncertainties describe the 
expected limits at 95% C.L. 


Solution 

Super-Kamiokande 

SNO 

SMA 

LMA 

1 S +4.2 

-'-•o -2.0 

7 8 +1-4 
-6.6 

9 -1-6.0 
^ -2.4 

15 +29 

-‘-a -13 


TABLE VI. What statistical test is best? The table compares the statistical power of three 
methods for analyzing data on the regeneration effect: a) moments of the zenith-angle distribu¬ 
tion, b) day-night asymmetry (An-d), and c) Kolmogorov Smirnov (K-S) test of the zenith-angle 
distribution. The number of sigmas listed in the table corresponds to the deviation of the best-fit 


MSW solutions described in Sec. from the undistorted zenith-angle distribution. The numerical 
results correspond to 30000 events for the SMA solution and 5000 events for the LMA solution. 




Moments 

An-d 

K-S test 

Detector 

Solution 


(o-) 

(o-) 

Super-Kamiokande 

SMA 

5 

4.4 

3.7 


LMA 

5.5 

9.2 

6.3 

SNO 

SMA 

6.5 

4.9 

4.4 


LMA 

10 

18.8 

12.4 
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TABLE VII. Moments of the Energy Spectrum. The moments in the absence of oscillations are 
To = 7.293 MeV and Cq = 3.391 MeV^ for Super-Kamiokande and Tq = 7.646 MeV and do = 3.032 
MeV^ for SNO. 


MSW Solution 

(T) 

(MeV) 

(T - To)/To 

(%) 

d2 

(MeV^) 

(%) 

Super-Kamiokande 

Day 

7.408 

1.58 

3.591 

5.88 

SMA 

Night 

7.403 

1.50 

3.574 

5.38 


Average 

7.405 

1.54 

3.582 

5.64 


Day 

7.275 

-0.25 

3.368 

-0.67 

LMA 

Night 

7.310 

0.23 

3.439 

1.42 


Average 

7.294 

0.008 

3.408 

0.48 


Day 

7.298 

0.06 

3.411 

0.58 

LOW 

Night 

7.290 

-0.04 

3.398 

0.20 


Average 

7.294 

0.008 

3.404 

0.38 

SNO 


Day 

7.869 

2.91 

3.161 

4.25 

SMA 

Night 

7.846 

2.62 

3.144 

3.68 


Average 

7.858 

2.76 

3.152 

3.96 


Day 

7.644 

-0.03 

3.032 

-0.03 

LMA 

Night 

7.719 

0.95 

3.113 

2.66 


Average 

7.687 

0.53 

3.084 

1.72 


Day 

7.682 

0.47 

3.069 

1.22 

LOW 

Night 

7.668 

0.28 

3.060 

0.90 


Average 

7.675 

0.37 

3.064 

1.06 
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TABLE VIII. Moments of the Energy Spectrum. The moments in the absence of oscillations are 
To = 8.057 MeV and Cq = 2.825 MeV^ for Super-Kamiokande and Tq = 8.178 MeV and ag = 2.348 
MeV^ for SNO. 


MSW Solution 

(T) 

(MeV) 

(T - To)/To 

(%) 

(MeV^) 

(%) 

Super-Kamiokande 

Day 

8.149 

1.14 

2.968 

5.07 

SMA 

Night 

8.143 

1.06 

2.954 

4.56 


Average 

8.146 

1.10 

2.962 

4.84 


Day 

8.046 

-0.14 

2.804 

-0.76 

LMA 

Night 

8.074 

0.22 

2.861 

1.25 


Average 

8.061 

0.05 

2.833 

0.28 


Day 

8.064 

0.08 

2.836 

0.40 

LOW 

Night 

8.058 

0.01 

2.827 

0.048 


Average 

8.061 

0.046 

2.832 

0.24 

SNO 


Day 

8.332 

1.88 

2.467 

5.05 

SMA 

Night 

8.315 

1.67 

2.452 

4.41 


Average 

8.323 

1.77 

2.459 

4.72 


Day 

8.177 

-0.043 

2.347 

-0.043 

LMA 

Night 

8.239 

2.63 

2.410 

2.63 


Average 

8.212 

1.84 

2.392 

1.84 


Day 

8.206 

0.34 

2.379 

1.28 

LOW 

Night 

8.196 

0.22 

2.370 

0.94 


Average 

8.201 

0.28 

2.374 

1.11 
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TABLE IX. Dependence of moments on assumed chemical composition. The table gives the 
relative shifts (in percent) of the first and second moments {A^i/fj,i) of the zenith-angle distribution 
in Super-Kamiokande and SNO as a function of the assumed chemical composition. The ratio ZjA 
in the core has been varied by ±0.5%, ±1%, ±2% (see second column) from the central value 
[(Z/yl)core = 0.465] adopted in the rest of the paper. The ratio in the mantle has been varied 
by —1%, —2% from the standard value of [(^/A)mantie = 0.496]. The last row for each detector 
corresponds to a simplified model with Zj A = 0.5 both in the mantle and in the core. The four 
columns for each detector correspond to the first and second moment in the SMA and LMA best-fit 
solutions respectively. 


Detector 

A(Z/A)/(Z/A) 

(%) 

SMA 

A(a)/ao 

(%) 

Aa^tal 

(%) 

LMA 

A(a)/ao 

(%) 

AaVcxg 

(%) 


-2 

1.03 

2.06 

3.22 

-0.338 


-1 

1.02 

2.05 

3.22 

-0.329 

Super-Kamiokande 

-0.5 

1.02 

2.04 

3.22 

-0.325 


0.0 

1.02 

2.03 

3.23 

-0.323 


±0.5 

1.02 

2.02 

3.23 

-0.316 

(core) 

±1 

1.01 

2.01 

3.23 

-0.311 


±2 

1.01 

1.99 

3.23 

-0.301 


-1 

1.01 

2.03 

3.19 

-0.336 

mantle 

-2 

1.00 

2.04 

3.16 

-0.358 

Z/A = 0.5 

0 

0.98 

1.89 

3.27 

-0.257 


-2 

1.05 

2.46 

5.27 

-2.23 


-1 

1.04 

2.44 

5.27 

-2.22 


-0.5 

1.04 

2.43 

5.27 

-2.21 


0.0 

1.04 

2.42 

5.28 

-2.20 

SNO 

±0.5 

1.04 

2.41 

5.28 

-2.20 
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+1 

1.03 

2.40 

5.28 

-2.19 


+2 

1.03 

2.38 

5.29 

-2.17 


-1 

1.01 

2.38 

5.23 

-2.15 

mantle 

-2 

0.98 

2.33 

5.18 

-2.10 

Z/A = 0.5 

0 

1.01 

2.28 

5.33 

-2.11 


TABLE X. Equatorial enhancement. The magnitude of the regeneration effect is compared for 
detectors located at their actual positions and at the equator. The size of the effect is represented 
by the number of standard deviations the first and second moments of the angular distribution 
differ from the undistorted exposure function (and, in parentheses, the day-night asymmetry in 
percent). The values of the neutrino oscillation parameters Am^ and sin^ 20 correspond to the 


best-fit SMA and LMA solutions (see Eqs. 1^,^ and 2a. For comparison, the best-ht SMA 
and LMA solutions produce 5 ct and ISu effects at the actual location of Super-Kamiokande and 
6.5(7 and 25 it at the actual location of SNO. 


Location 

Solution 

Super-Kamiokande 

SNO 

Equator 

SMA 

LMA 

«<r (3.0 %) 

14.5,7 (7.6 ^^(*3 %) 

9ct(4.1 +li%) 

29a (15 tit %) 
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FIG. 1. Schematic view of detector’s location and sun’s direction. The zenith is defined as the 
line from the center of the Earth through the center of the detector. The zenith angle, a, and the 
latitude of the detector’s location, (f), are also shown in the figure. 


FIG. 2. Density profiles for six different models of the earth. The models are: I) PREM |37|, 


2) Stacey’s model |Q, 3) model HA Q, 4) model HBl Q, 5) model B497 |Q, and 6) model 
B1 [42|. The radius is given in kilometers from the center of the earth. All models are spherically 
symmetric. 


FIG. 3. Allowed MSW solutions with regeneration. The allowed regions are shown for the 
neutrino oscillation parameters Am^ and sin^ 20. The G.L. for the outer regions is 99% and the 
C. L for the inner regions is 99% (only applies to the LMA and SMA solutions). The data used 
here are from the chlorine |j^,Kamiokande Q,GALLEX Q, and SAGE |Q] experiments. The solar 
model used is the best standard model of Bahcall and Pinsonneault (1995) with helium and heavy 
element diffusion Q. The points where has a local minimum are indicated by a circle. 


FIG. 4. Survival probabilities for MSW solutions. The figure presents the survival probabilities 
for a t'e created in the sun to remain a upon arrival at the Super-Kamiokande detector. The 


best-fit MSW solutions including regeneration in the earth are described in Sec. IV. The full line 
refers to the average survival probabilities computed taking into account regeneration in the earth 
and the dotted line refers to calculations for the day-time that do not include regeneration. The 
dashed line includes regeneration at night. There are only slight differences between the computed 
regeneration probabilities for the detectors located at the positions of Super-Kamiokande, SNO 
and the Gran Sasso Underground Laboratory. 
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FIG. 5. Super-Kamiokande and SNO zenith-angle distributions. The figure shows the expected 
zenith-angle distribution of neutrino events during one calendar year in the Super-Kamiokande and 
the SNO detectors. The angle, a, represents the angular separation between the direction to the 
sun and the direction of the local zenith (see Fig. ||). The two left panels display the zenith-angle 
exposure functions, the undistorted angular distributions in the absence of oscillations. The ex¬ 
posure functions are determined by the location of the two detectors at, respectively, Kamioka, 
Japan, and Sudbury, Canada. The distorted zenith-angle distributions due to the regeneration 
effect in the Earth are shown in the two right panels; the neutrino solutions are indicated by: SMA 
(solid line) and LMA (dotted line). 

FIG. 6. Gran Sasso zenith-angle distributions. The figure shows the expected zenith-angle 
distribution of events during one calendar year in detectors located at the Gran Sasso Laboratory in 
Italy: ICARUS, BOREXINO, HERON and HELLAZ. The upper left panel shows the zenith-angle 
exposure function, which does not depend on detector characteristics. The three additional panels 
display the distorted zenith-angle distributions due to the regeneration effect in the Earth; the 
solutions are indicated by: SMA (solid line), LMA (dotted line) and LOW (dashed line). 

EIG. 7. Contours of constant relative shift (in percent) of the average zenith angle, ((a)—Q;o)/aO) 
due to t'e regeneration in the earth as a function of the neutrino oscillation parameters, sin^ 29 
and Am^. Here uq = 90^ is the average angle of the undistorted angular distribution with no 
oscillations. The shaded regions in the panels for Super-Kamiokande and SNO are allowed by the 
latest solar neutrino data at 95% C.L. and represent the SMA and LMA solutions. In the lower two 
panels (BOREXINO and HERON/HELLAZ) the three shaded regions are allowed at 99% C.L., 
the low-mass region representing the LOW solution (see text for details). The black circle within 
each allowed region represents the point which corresponds (locally) to the best-fit to the data. 
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FIG. 8. Contours of constant relative shift (in percent) of the dispersion of the zenith angle, 
(d^— cjq)/(Tq, due to Vf. regeneration in the earth as a function of the neutrino oscillation parameters, 
sin^ 20 and Am^. The values of ctq are given in the text for each of the experiments. The definition 
of the shaded regions is the same as in Fig. 


FIG. 9. How many sigmas? The figure shows the sensitivity of Super-Kamiokande, SNO, 
ICARUS, BOREXINO and HERON/HELLAZ to the regeneration effect. Iso-sigma contours, sta¬ 
tistical errors only, delineate the fractional percentage shifts of the first two moments of the angular 
distribution of events for an assumed 30000 observed events. Eor all but the ICARUS experiment, 
the best-ht MSW solutions are indicated by black circles (SMA), squares (LMA), and triangles 


(LOW); the best-fit solutions are presented in Sec. W. The error bars on the predicted moments 
correspond to Am^ and sin^ 20 within allowed solution space at 95% C.L. (for Super-Kamiokande, 
SNO, and ICARUS) or 99% C.L. (BOREXINO and HERON/HELLAZ). Eor ICARUS, we have 


indicated the best-fit solutions by a transparent circle, square, or triangle. The best-fit SMA and 
LOW solutions for ICARUS and the LOW solution for SNO are all three close together at about 
3a from the no oscillation solution. In order to avoid too much crowding in the figure, we have not 
shown the theoretical uncertainties for ICARUS. 


FIG. 10. Contours of constant day-night asymmetry, An-d (see Eq. (^)), in Super-Kamiokande, 
SNO, ICARUS, BOREXINO and HERON/HELLAZ. The shaded regions are the same as in Fig. 0. 


FIG. 11. Relative distortion for Super-Kamiokande. The figure shows the fractional distortion, 
[/(a) — Y{a)]/Y{a), of the zenith-angle distribution for the best-estimate SMA and LMA MSW 
solutions. 

FIG. 12. The zenith-angle distribution for equatorial detectors with the characteristics of Su¬ 
per-Kamiokande, SNO, ICARUS, BOREXINO and HERON/HELLAZ. Notation is the same as in 
Figs. ^ and 
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FIG. 13. How many sigmas at the equator? The figure shows the sensitivity to the regenera¬ 
tion effect of equatorial detectors with the characteristics of Super-Kamiokande, SNO, ICARUS, 
BOREXINO and HERON/HELLAZ. Notation is the same as in Fig. p. 
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